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ABSTRACT 

Using the NIRSPEC spectrograph at Keck II, we have obtained H and iGband echelle 
spectra for a young (~ 10 — 15 Myr), luminous (My ~ —13.2) super-star cluster in the 
nearby spiral galaxy NGC 6946. From spectral synthesis and equivalent width mea- 
surements we obtain for the first time accurate abundances and abundance patterns 
in an extragalactic super-star cluster. We find [Fe/H]= —0.45 ± 0.08 dex, an average 
a- enhancement of » +0.22 ± 0.1 dex, and a relatively low 12 C/ 13 C » 8 ± 2 isotopic 
ratio. We also measure a velocity dispersion of «9.1 km/s, in agreement with previ- 
ous estimates. We conclude that integrated high-dispersion spectroscopy of massive 
star clusters is a promising alternative to other methods for abundance analysis in 
extragalactic young stellar populations. 

Key words: Galaxies: individual (NGC 6946), star clusters, abundances — infrared: 
galaxies — techniques: spectroscopic 



1 INTRODUCTION 

Star clusters have a long history as important tools for ad- 
dressing a wide range of questions in astronomy. With few 
exceptions, they are "simple stellar populations" (SSPs), i.e. 
they are composed of stars born in a single burst and shar- 
ing the same chemical composition (at least to first order) 
and age. They have played an i mportant role as test labs for 
models of stellar evolution ( e.g. iMaeder fc Mermilliodlll98lt 
iRenzini fc Fusi Peccil ll98Sfe IChiosi et all Il992l) . As confi- 
dence has grown in our ability to model their integrated 
properties, so has their importance as tracers of stellar pop- 
ulations in galaxies that are too distant for individual stars 
to be resolved. 

The latter point is perhaps best illustrated by consider- 
ing, for the moment, the old globular clusters (GCs) which 
are ubiquitous in all major galaxies. GCs typically contain 
large numbers of stars (10 4 — 10 6 ), so most phases of stellar 
evolution are well sampled and stochastic effects therefore 
minimized. Although there are still unsolved problems (e.g. 
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concerning horizontal branch morphology), SSP models can 
now provide a reasonably realistic account of integrated GC 
observables such as broad-band colours and absorption line 
strengths as a function of age and metallicity. Thanks to 
the availability of efficient spectrographs on 8-10 m class 
telescopes, spectroscopy of GCs in galaxies well beyond the 
Local Group is now routinely feasible, and has been utilized 
in many studies to probe the stellar populations in early- 
type galaxies and const rain their ages and metallicities (e.g. 
iBrodie fc Straderl2006l . and references therein) . Much of this 
work has relied on measurements of absorption line features 
at relatively low (~ 10 A) spectral resolution. In the op- 
tical, the Lick/IDS syst em of absorption line indices has 
found wide-spread use llBurstein et al.l fl984t iTraeer et alJ 
1998) , while indic e s for t he ne ar-infrared have been defined 
bv iBica fc Alloinl dl987l) and IVazdekis et"ai1 (120031) . Ages 
are typically estimated from Balmer line strengths (mainly 
H<5, H7, and H/3), with other indices being more sensitive 
to metallicity. However, detailed abundances of individual 
elements typically cannot be reliably measured at this reso- 
lution. In addition, the Lick/IDS system is not well tailored 
for studying young stellar populations, partly due to lim- 
itations in the empirical libraries used in the construction 
of SSP models, but also because the index definitions are 
designed primarily for studies of old stellar populations. 

Until now, the main source of information about the 
chemical composition of young stellar populations beyond 
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the Local Group has been measurements of emission lines in 
Hn regions. However, these provide access to only a limited 
set of elements (mainly O, N, S) and often have to rely 
on empirical calibrations of line strengths vs. metallicity. 
The auroral lines (e.g. [Om] 4363A) are usually too faint 
to be measured directly, t hus prohibiting a determination of 
the nebular temperature (Stasiriska 2001). Furthermore, Hll 
regions only offer a snapshot of the present-day chemistry 
and do not provide any information about the past history. 

Spectroscopy of (massive) star clusters has the poten- 
tial to provide information on the entire star formation his- 
tories of galaxies. For masses in the 10 s M© - 10 6 Mq range, 
star clusters typically have velocity dispersions of about 5- 
10 km/s, allowing studies of their integrated properties at 
spectral resolutions up to A/ A A ft! 30,000 or more. This is 
sufficient to constrain individual element abundances. How- 
ever, extending this type of analysis to young clusters comes 
with its own set of difficulties. Compared to studies of old 
GCs, some of the main challenges are: 1) models for the 
massive stars found in clusters with ages of a few x 10 7 years 
are much les s certain than those of t he low-mass stars found 
in old GCs iMassev fc Olser] 120031) 2) the optical spectra 
are generally dominated by the relatively featureless spec- 
tra of hot stars, diluting absorption features and making 
them harder to measure, 3) for all but the most massive 
clusters, the integrated spectra may be dominated by a few 
massive stars, causing unpre dictable stochastic fluctuat ions 
in integrated properties fe.g. lLancon fc Mouhcinell2000l) . 

In this letter we aim to take a first step towards deriv- 
ing abundances for young extragalactic star clusters from 
their integrated light, by modelling the near-infrared spec- 
trum of a luminous (My = —13.2) young (~ 10 Myr) star 
cluster in the nearby spiral galaxy NGC 6946. This object 
was fi rst identified as a star cluster by lLarsen fc Richtlerl 
Jl999tl and was labelled NGC6946-1447 by them, b ut is lo- 
cated within a larger stellar complex first noted by iHoded 
lll967l) . Both UBVI broad-band colours and spectroscopic 
observations of Balmer and He I absorption li nes yield con- 
sistent age estimates in th e range 10-15 Myr jLarsen et alJ 
l200ll : lEfremov et aljl2002l) . From the luminosity and esti- 
mated age of the cluster, SSP models yielded a mass of 
~ 0.8x 10 6 M for a Salpeter-like IMF extending down to 0.1 
M , or ~ 0.55 x 10 6 M w if the IMF is log- normal below 0.4 
Mq (lLarsen et al.ll200ll) . A dynamical mass estimate, based 
on high-dispersion optical spectroscopy and HST imaging, 
yielded a somewhat higher mass of ~ 1.7 x 10 6 M . Such 
massive clusters are sometimes called "super star c lusters" 
(SSCs). Using isochrones from iGirardi et al] feOOdL wc es- 
timate that the cluster contains about 130 red supergiants. 
A detailed descripti on of the clus t er and surro unding stellar 
complex is given in lLarsen et al] i200 li . I2002T) . 

No line emission is observed from the cluster itself but 
lEfremov et al] i2002T) estimated an oxygen abundance of 
12 + \og(0/H) = 8.95 ± 0.2 from long- slit spectroscopy of 
nearby Hll regions. iBellev fc Rovl (1992) measured O abun- 
dances for Hll regions distributed throughout the disk of 
NGC 6946 (using narrow-band imaging of O, N and H emis- 
sion lines) and derived an oxygen abundance gradient of 
Alog(0/H)/AR = -0.089 ± 0.003dexkpc _1 and a central 
value of 12 + log(0/i/) = 9.37. At the projected galactocen- 
tric distance of NGC6946-1447 (4 8 kpc), this corre s ponds 
to 12 + \ag{0/H) = 8.94 ± 0.01. iKobulnickv et al] Jl998l) 



quote an O abundance of 12 + log(0/_ff) = 9.13 at a radius 
of 3 kpc, or 12 + \og{0/H) = 8.97 at 4.8 kpc if w e adopt 
the abundance gradient from IBellev fc Rovl Jl992» . Thus, 
all available measurements consistently give 12 + \og(0/H) 
between 8.94 and 8.97, although the above O abundances 
are all based on measurements of the collisionally excited O 
lines vs. Balmer line ratios and may therefore be subject to 
systematic uncertainties at the 0.1-0.2 dex level. 



2 OBSERVATIONS AND DATA REDUCTION 

H and K-band high-resolution spectra were acquired on 
13 July 2002 using the infrared spectrograph NIRSPEC 
iMcLean et al]ll998h mounted at the Nasmyth focus of the 
Keck II telescope. The high resolution echelle mode, with a 
slit width of 0V43 (3 pixels) and a length of 12" was used, 
providing a spectral resolution of A/AA = 25,000. The ex- 
posure times were 64 min and 48 min in the H and K bands, 
respectively, yielding a S/N of about 26 and 36 per pixel in 
the dispersion direction (the K magnitude of the cluster is 
about 13.0). The observations were obtained in pairs of ex- 
posures with a duration of 240 s each, nodded a few arcsec 
along the slit to allow reliable sky subtraction without any 
additional overhead for separate sky e xposures. The NIR - 
SPEC data has previously been used bv lLarsen et al] (120041) 
to derive a line-of-sight velocity dispersion of 8.8 ± 1.5 km/s 
for the cluster, i n good agreement w ith the value of 10.0±2.7 
km/s derived bv lLarsen et al.H200l]) (based on Keck/HIRES 
spectroscopy). We referlio TLarseTietaTI i2004f) for details on 
the data reduction. 



3 SPECTRAL ANALYSIS 

Near-IR spectroscopy is a powerful tool to obtain accurate 
abundances of key metals in cool stars (T e ff < 5000 K). Sev- 
eral atomic and molecular lines are strong and not affected 
by severe blending, making them powerful abundance trac- 
ers not only in stars but also in more distant stellar clus- 
ters and galaxies and for a wide range of metallicities and 
ages dOriglia et al.ll997l : lOliva fc Origliell99gt lOriglia et all 
2004). However, to properly account for line blending, abun- 
dance analysis from integrated spectra generally still re- 
quires full spectral synthesis techniques and not just equiv- 
alent width measurements of individual lines. Population 
synthesis may also be required to define the dominant con- 
tribution to the stellar luminosity. 

The near IR stellar continuum of young stellar clusters 
and starburst g alaxies is almost enti rely due to luminous 
red supergiants JOriglia fc Olivall2000l) and usually it domi- 
nates over nebular and dust emission. Based on the Girardi 
et al. isochrones, stars hotter than 10,000 K contribute only 
~ 5% of the H- and K-hanA flux at 15 Myrs. This repre- 
sents a major, conceptual simplification in population and 
spectral synthesis techniques, making the interpretation of 
the integrated spectra much easier. The spectra can be mod- 
elled with an equivalent, average star, whose stellar param- 
eters (temperature T e ff, gravity log<? and microturbulence 
velocity £) mainly depend on the stellar age and metallicity. 
Both observations and evolutionary models (see e.g. iKellerl 
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Figure 1. Near-IR spectra of the SSC of NGC 6946. Observed spectra: histograms; synthetic stellar best-fit solution: solid lines. A few 
atomic and molecular features of interest are also marked. 



Il999l: lOrielia et al1ll999l : iMassev fc rilsenll200l and refer- 
ences therein) suggest that red supergiants of ages between 
~6 and 100 Myr and metallicities between 1/10 and Solar 
have low gravities (log g<1.0), low temperatures (< 4000 K) 
and relatively high microturbulence velocity (£ >3 km/s). 

At the NIRSPEC resolution of R=25,000, several sin- 
gle roto-vibrational OH lines and CO bandheads can be 
measured and used to derive accurate oxygen and carbon 
abundances. Although our NIRSPEC setup was optimised 
for velocity dispersion measurements rather than abundance 
analysis, abundances of other metals can be derived from the 
atomic lines of Fe I, Mg I, Si I, Ti I, Ca I and Al I. 

A grid of synthetic spectra of red supergiant stars for 
different input atmospheric para meters and abundances 
were computed, using an updated dOriglia. Rich fc Castrol 
2002; Oridia ct al. 2003) version of the code described 
in lOriglia. Moorwood fc Olival 1(1993). Briefly, the code 
uses the LTE approximation a nd is based on molecular 
blanke ted model atmospheres of I.Tohnson. Bernat, fc KruprJ 
(1980) at temperatures <4000 K and the ATLAS9 models 
for temperatures above 4000 K. Recently, the NextGen 



model atmospheres ijHauschildt et al"1ll999l) have been also 
implemented within the code and tested. Compared with 
the older models, the differences in the resulting abun- 
dances are only min or (well within a few hundredths dex; 
iRich fe Origlial2005tl . This is not surprising, since the major 
source of opacity in the near IR spectra of cool stars is H~ 
with a minimum around 1.6/im and small differences in the 
temperature structure of different model atmospheres have 
a minor impact on the overall abundance determination. 
The code also includes several thousands of near IR atomic 
lines and molecular roto-vibrational transitions due to CO, 
OH and CN. Three main compilations of atomic oscilla- 
tor strengths are used, namely the Kurucz's database (c.f. 
http:/ / cfa-www. harward. e du/ amdata/ ampdata/ kurucz23/sekur. html I , 
and those published bv i Biemont fe Grevessd lll97Stl and 
iMelendez fc Barbuvl l|l999E 

The code provides full spectral synthesis over the 1-2.5 
/im range and abundance estimates are mainly obtained by 
best-fitting the full observed spectrum and by measuring the 
equivalent widths of a few selected features (Fig. , dom- 
inated by a specific chemical element, as a further cross- 



4 Larsen et al. 

Table 1. Adopted stellar atmosphere parameters and abundance estimates for the SSC in NGC 6946. 
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a [a/Fe] is the average [< Ca, Si, Mg, Ti > /Fe] abundance ratio. 
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Figure 2. Average probability of a random realization of our 
best-fitting solution and the test models with varying temperature 
by AT eff of ±200K (left panels), microturbulence by A§ of :p 1.0 
km s -1 (right panels), and gravity by Alog g of ±0.5 dex (middle 
panels), with respect to the best-fitting solution (see Sect.[3J for 
the SSC. 

check. The equivalent widths were measured by perform- 
ing a Gaussian fit with a equal to the measured stellar 
velocity dispersion; typical values range between 100 and 
500 mA with a conservative error of ±20 mA to also ac- 
count for a ±2% uncertainty in the continuum positioning. 
By best-fitting the full observed IR spectrum and by mea- 
suring the equivalent widths of selected lines, we obtained 
the following stellar parameters and abundance patterns: 
T cff =4000, log 5 =0.5, £=3, [Fe/H]=-0.45; [O/Fe]=±0.28, 
[<Si,Mg,Ca,Ti>/Fe]=±0.22; [Al/Fe]=±0.25; [C/Fe]=-0.25. 
We also measure 12 C/ 13 C » 8 ± 2. 

Table lists the derived abundances and their asso- 
ciated random erro rs at 90% confidence. Re ference Solar 
abundances are from lGrevesse fe Sauvall jl998l) . In addition, 
we also measure a heliocentric radial velocity of < v r >= 
±141 ±2 km/s and velocity dispersion (corrected for instru- 
mental broadening) of « 9.1 ±1 km/s. These numbers agree 
reassuringly well with previo us estimates iEfremov et alJ 
l2002l:lLarsen et alJl200lLl2004h . 

3.1 Error budget 

Synthetic spectra with lower element abundances are sys- 
tematically shallower (have weaker features) than the best- 
fit solution, while the opposite occurs when higher abun- 



dances are adopted. In order to check the statistical signifi- 
cance of our best-fit solution, as a function of merit we adopt 
the difference between the model and the observed spectrum 
(hereafter S). This parameter is more effective for quantify- 
ing systematic discrepancies than the classical \ 2 test, which 
is equally sensitive to random and systematic scatters (see 
lOriglia et aT] l2003. for more details). 

Since 5 is expected to follow a Gaussian distribution, 
we compute S and the corresponding standard deviation 
(a) for the best-fit solution and 6 test models with stel- 
lar parameters varying by AT c ff±200 K, Alog g=±0.5 and 
A£=pl.O km s~ with respect to the best-fit, and abundances 
varying accordingly by ~ ±0.2 dex, in order to still repro- 
duce the depth of the observed features. We then extract 
10000 random subsamples from each test model (assuming 
a Gaussian distribution) and we compute the probability P 
that a random realization of the data-points around a test 
model display a S that is compatible with the best-fit model. 
P ~ 1 indicates that the model is a good representation of 
the observed spectrum. 

The left panel of Fig. |5] shows the results for the ob- 
served H band spectrum of the SSC in NGC 6946. It can be 
easily appreciated that the best-fit solution provides in all 
cases a clear maximum in P (>99%) with respect to the test 
models, which are statistical significant only at > 1.5a level. 
We also computed test models with the same stellar parame- 
ters as the best-fit solution and varying only the abundances. 
Models with ±0.1 dex with respect to the best-fit solution 
are significant at 1 — 1.5<7 level, while models with ±0.2 dex 
are only marginally acceptable at a > 3a level. 

Hence, as a conservative estimate of the systematic error 
in the derived best-fit abundances, due to the residual uncer- 
tainty in the adopted stellar parameters, one can assume a 
value of « ±0.1 dex. Moreover, since the stellar features un- 
der consideration show a similar trend with variation in the 
stellar parameters, although with different sensitivity, rela- 
tive abundances are less dependent on stellar parameter as- 
sumptions, reducing the systematic uncertainty to <0.1 dex. 



4 DISCUSSION AND CONCLUSIONS 

Based on near-infrared H and if-band spectroscopy, we have 
carried out a detailed abundance analysis of a young massive 
star cluster in the nearby spiral galaxy NGC 6946. In this 
exploratory work, we have derived abundances of several 
individual key elements, including Fe, Ca, Si, Mg and Ti. 
We find a sub-solar Fe abundance ([Fe/H] = -0.45 ± 0.08), 
while the a-element to Fe abundance ratios are all enhanced 
with respect to the Solar values with a mean [a/Fe] = 
0.22±0.11. The O abundance derived from the cluster spec- 
trum is [0/H]= -0.17 ± 0.09, about 0.3 dex lower than the 
value based on Hn regions at the same galactocentric dis- 
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tance. We find a 12 C/ 13 C ratio of « 8 ± 2, similar to or 
slightly lower than typi cally observed in red supergiants in 
Gala ctic open clusters ()Lucklll994t lOonzalez fc Wallersteinl 
2000) and in the (more metal-poor) cluster NGC 330 in the 
SM C (iGonzalez & , Wallersteiril999T) . Standard stellar mod- 
els JSchaller et alJll992ft predict a surface 12 C/ 13 C ratio of 
~ 17 for a 15 M Q , ~ 12 Myr-old (Solar metallicity) red su- 
pergiant, about a factor of two higher than the value derived 
here. 

Super-solar [a/Fe] ratios are usually interpreted as 
signatures of rapid, bursty star formation, with the gas 
mainly enriched by Type II supernovae w ith short-lived, 
massive progenitor stars jMcWilliarr]ll997ft . Stellar popu- 
lations formed over timescales of several Gyrs or more are 
expected to show solar-like a/Fe-element abundance ratios, 
as observed in the Milky Way thin disk. 

In this context, it is worth noting that the complex 
hosting NGC6946-1447 may qualify as a 'localized starburst' 
jEfremovll2004h . Timing is a critical issue, however: in order 
to produce super-solar [a/Fe] ratios, the starburst must have 
preceded the formation of the cluster itself by several Myrs. 
The delay must have been long enough for Type II SNe to 
produce significant amounts of Q-elements, and enough time 
must then have elapsed to allow mixing of the pre-existing 
gas with a-enhanced ejecta before the cluster formed. Inter- 
estingly, reconstruction of the field star colour-magnitude 
diagram has provided some evidence for star formation in 
the complex at least 10-15 My r prior to the formation of 
the cluster jLarsen et ail fecx^. Alternatively, the current 
global SFR in NGC 6946 may be sufficiently elevated above 
the past average to cause a general enrichment of the ISM 
with a-elements. Addressing these questions more quanti- 
tatively would require a detailed modeling of the chemical 
enrichment history and a knowledge of the past star forma- 
tion history which is currently not available. 

This study represents one of the first cases of a detailed 
abundance analysis of a star-forming galaxy beyond the Lo- 
cal Group. While many uncertainties remain, we suggest 
that observations of extragalactic young star clusters hold 
great potential for constraining the chemical enrichment his- 
tories of their parent galaxies. 
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